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CN ■ Abstract 

> _ 

In the present paper we study the Faddeev-Popov path integral quan- 
tization of electrodynamics in an inhomogenious dielectric medium. We 
quantize all polarizations of the photons and introduce the corresponding 
ghost fields. Using the heat kernel technique, we express the heat kernel 
coefficients in termini of the dielectricity e(x) and calculate the ultra violet 
divergent terms in the effective action. No cancellation between ghosts and 
"non-physical" degrees of freedom of the photon is observed. 



On 



-F; ■ 1 Introduction 



The Casimir effect describes the forces resulting from the vacuum fhictua- 
$_i ■ tions (ground state energy) of the electromagnetic field in simple situations 

realized by conducting surfaces. These forces can be viewed as retarded 
Van der Waals forces between the atoms constituting the surfaces (and 
the bodies behind). As a generalisation of this picture one can consider 
some medium. It can be characterised either by atoms at positions Xi 
with their individual polarisabilities on or by a macroscopic permittivity 
e(x) and permeability [x(x). Again, we can calculate the resulting poten- 
tial of the VanderWaals forces or the vacuum energy Eo[e(x) , fj,(x)] of the 
electromagnetic field in a background given by e(x) resp. n(x). Taking 
into account that real permittivity resp. permeability are functions of the 
photon frequency we arrive at the problem to calculate Eo[e(x, uj), /i(x, u)]. 
The dependence on u> has as a physical background, besides others the ob- 
servation that any medium becomes transparent for u sufficiently high (we 
do not consider inelastic effects here). Therfore e, \x — ► 1 for uj — > oo should 
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serve as a natural ultraviolet regularisation. This is widely believed, but 
not shown in a rigorous way yet. 

The problem of the calculation of Eo[e(x), fJ,(x)], i.e. without frequency 
dependence, may be well posed independently. A physical justification 
could be that the essential contribution results after a proper renormal- 
isation from quite low frequencies uj, where e and \i can be viewed as 
approximately independent on u. In that case we don't have a natural 
regularisation and have to proceed like in the general situation with sharp 
boundary conditions or a general background field. For technical reasons 
we use the zeta-functional regularisation. Then the first step is to calculate 
the divergent contributions (the proper technical tool being the heat ker- 
nel expansion), the second is to formulate a model for the interpretation 
of the renormalisation (this is to be able to reinterpret the subtraction of 
the divergencies as a renormalization of classical quantities like volume, 
surface tension etc. as discussed in Q or like mass and coupling constant 
of the background field as discussed in P]) and, finally, in a third step to 
calculate the renormalised groundstate energy Eq . In the present paper we 
carry out the first step and discuss the second to some extend. 

The forces resulting from the electromagnetic vacuum fluctuations in 
polarisable media have been given much attention to. The common fea- 
tures of these investigations are sharp boundaries separating regions of 
different values of e{x) and simple geometries (planes, cylinders and so 
on). For instance, there have been interesting calculations (mainly with 
respect to the sign of the force) for a dielectric sphere ||. Also, much at- 
tention had been spent on a possible explanation of sonoluminescence as a 
dynamical Casimir effect, especially in a series by Schwinger [Q]. Recently, 
the bulk and surface energy contributions had been discussed [||, ||]. 

However, with respect to the renormalisation it is difficult to deal with 
sharp boundaries resp. non smooth background fields. It is known that 
additional contributions to the heat kernel expansion occur and therefor 
additional counter terms result for which a general theory is still missing. 
Therefor we restrict ourselves in the present paper to e(x) which are smooth 
functions on x. 

There is still another problem we have to pay attention. In the common 
understanding the quantisation of QED in media is done in the Coulomb 
gauge, i.e., the two 'physical' polarisations of the photon are quantised. 
Also there are known procedures where all polarisations of the photon 
are quantised and the gauge invariance (in the presence of boundaries) 
is restored by ghosts which have to fulfil boundary conditions too (one 
of the first is 0], later on it had been discussed in Q). In most cases 
their contributions cancel that resulting from the 'unphysicak photons, but 
counter examples are known (e.g. for QED in curved space time, ||). In 
the framework of quantum optics the canonical quantization of photons was 
considered in [|ll| without, however, analysing the ghost contributions. An 
alternative approach for quantization in covariant gauge without ghosts, 
however restricted to sharp boundaries, had been developed in [10]. 
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In the present paper we analyse the problem of QED with a position 
dependent permittivity e{x) from the point of view of general quantum 
gauge theory in an external field. We analyse the canonical path integral 
measure and corresponding configuration space measure. A gauge fixing 
term is introduced together with the ghost action. Next we analyse the 
ultra violet structure of the theory by means of the heat kernel expansion. 
No cancellation between ghosts and photon modes is obtained. 

Our paper is organized as follows. In the next section the quantization 
of the theory is considered and the path integral is derived. In sec. 3 we use 
the heat kernel expansion to evaluate ultra violet divergencies. Concluding 
remarks are given in sec 4. An Appendix contains an alternative calculation 
to check up the results of sec. 3. 



2 Canonical quantization and gauge choice 

Consider the action for the electromagnetic field in a dielectric media with 
permittivity e{x): 

S = J *x\{e{x)& - &) (1) 

To avoid technical complexities we put the permeability fi = 1 and suppose 
that e depends on spatial coordinates only. 

Let us rewrite the action (|T|) in the canonical first order form: 

Si = / rfxiFdoAi + A diP* - -^P'P 1 - \b 2 ) (2) 
J e(x) 2 

Here An is the vector potential. P l = —e{x)Ei is the momentum conjugate 
to A{. Canonical Poisson brackets are 

{A(x,t),P J '(y,t)} = ^(x-y) (3) 



The same brackets were obtained in Aq plays the role of a Lagrange 
multiplier generating the Gauss law constraint, which in turn generates 
gauge transformations. According to the general method 12] of quantiza- 
tion of gauge theories we can write down the path integral 

Z = J DAiDAoDPijppSixiAVeMiSi), (4) 

where x(Ai) is a gauge fixing condition, Jpp is the Faddeev-Popov de- 
terminant, Jpp = det{x(A), djP 3 }. Now we can perform the integration 
over the momenta PK It produces the factor ]J X \f e ( x ) A which should be 
absorbed in the path integral measure DA{. We arrive at the following 
expression: 

Z = J DAiDAoJFp5(x(A))exp(iS), M = yfeU (5) 
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Our A variables coincide with the q' variables of Glauber and Lewenstein 
Q. Note that the measure in (||) differs from the naive one Y\DA^. We 
can use the Faddeev-Popov trick to transform the path integral (||) to 
whatever gauge condition we prefer, introduce a gauge fixing term and 
ghost fields. There is nothing specific in this respect in the present model. 
All steps repeat those of a standard text book [12]. The result is 

Z = J DAiDA Dc Dcexp ji J d A x[^(2e(x)(d e- 1/2 Ai - diA ) 2 

-(cV" 1/2 i fe - 46- 1 / 2 !,) 2 ) + C gf + C ghost }} (6) 

where C g f and Cghost are gauge fixing term and ghost action respectively. 
As usual, we can bring the action in (0) to the form / A fl L^ tu A u , where 
Lfiv is a second order differential operator. In calculating the effective 
action and the heat kernel expansion it is much more convenient to deal 
with operators of Laplace type, i.e. operators with scalar leading symbol. 
There is a unique gauge choice which splits the in a direct sum of 
operators of Laplace type. This choice is 

C gf = ~(e- 1 d l e 1 / 2 A l -ed A ) 2 (7) 

Cghost = -c{-t~ l ditdi + edl)c (8) 
The action for the electromagnetic field A then takes the form 

l - J d^diAof - e 2 (d A ) 2 + (doA,) 2 (9) 

+Aie~ 1/2 (d 2 5 ik - eid k + die k - eie k )e~ 1/2 A k ], a = d^lne 

Note, that the mixing between Aq and Ai is removed completely. 

The total action with gauge fixing and ghost term is invariant under 
the BRST transformations with the parameter a(x): 

5Aq = done 

5Ai = e 1/2 diac 

5c = (10) 

5c = (-e~ 1 d i €^ 2 A i + ed A )a 

which are given here to complete the picture. 

3 Effective action and heat kernel expan- 
sion 

Now we are able to integrate over Aq, A and the ghosts. The resulting 
path integral reads after Wick rotation to the Euclidean domain: 

Z = Z[Aq]Z[A]Z[c,c], (11) 
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where the separate contributions are of the form: 
Z[A ] = dzt-^i-diedi-^dl) 

Z[A] = det- 1 / 2 f—dlSij-d^Sij-Gidj + Gjdi-Mij) (12) 

Z[c,c] = det{-e- l diedi- edl) 

We introduced the notations: 

e- 1 

G i = ~ Mij = -(dj - eiej) (,j Oji , (13) 

For the functional determinants we use the integral representation 

r°° dt 

logdet(L)= / —K(L; t) (14) 
Jo t 

where the heat kernel K(L; t) for a second order elliptic operator L is 

K(L;t) = Trexp(-tL) (15) 

The ultraviolet behavior of functional determinants is given by the 
asymptotic expansion of the heat kernel (15) as t — > +0. Since all the 
operators are of Laplace type, we can use the general theory fP3fl . Each of 
the operators has the structure 

L = -{ g rd^d v + a a d tr + b) (16) 

where g^ v plays the role of a metric. a u and b are local sections of endo- 
morphism End(y) of certain vector bundle. By introducing a connection 
•jjp, in the vector bundle V, one can bring L to the form: 

L = -(g^V t y u + E) (17) 

where V is a sum of the Riemannian covariant derivative with respect to 
the metric g and the connection u. The explicit form of u and E is 

us = \g v 5{a» + g^ a ) 

E = b-g^id^ + u^-uj^) (18) 

As usual, r denotes the Christoffel connection. 

Given the geometric quantities g, uj and E, we are able to calculate the 
coefficients a n of the asymptotic expansion 

oo 

Tr(/ exp(-tL)) = t~ 2 ]T t n a n (f, L) (19) 

n=0 

for a function /. The coefficients a n (f,L) contain information on the 
asymptotics of the heat kernel diagonal < x\ exp(— tL)\x >. The analytical 
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expressions for the first coefficients are known [JDJ: 
ai = j^trv J 'dW /2 f{E+ T ~) 

a 2 = j^y^v J dW /2 3^o/( 60£J ;/ + 60^ + 180^ + 30^^ 

+ 12T ; / + 5T 2 -2p 2 + 2i? 2 ) (20) 

Here -R, p and r are Riemann tensor, Ricci tensor and scalar curvature 
of the metric g respectively. Semicolon denotes covariant differentiation, 
E.p = Vu-E- All indices are lowered and raised with the metric tensor, try 
is the bundle (matrix) trace, VI is the field strength of the connection u: 

= d^uj v - dvUn + uj^uj v - (21) 

The three coefficients ( p0| ) are enough to describe the one-loop ultra 
violet divergencies in four dimensional quantum field theory in an infinite 
space-time. 

Now our problem is reduced to the calculation of the geometric quan- 
tities appearing in (^). For the ghost operator we have 

9ij = 9oo = ^(x) 

1 oo — 2 e % 0i ~ 2 
3 

ujo = 0, ui = -a 

3 _3_ 

4 16 

R'jki = (22) 

R°i0j = — 2 6 * J 7 e * e -? 
1 1 

Poo = ^(eii - 2 e i e i) 
1 

t — an 2^*^* 
For the operator acting on Aq the relevant quantities are: 
9ij = 500 = ^ 2 {x) 

1 00 — e « 1 0i — e « 

5 

ojq = 0, cjj = -ej 
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R jkl = -^(—ejidik + ejkdu + eu5kj — eikSij) (23) 
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+-{e p e p (5ji5 k i - 5j k Su) + e k ejSu - e k ei5ji - eiej5 k i + ei^jk) 

R°iOj = — e-ij + -jZie-ihj 

3 1 1 _3 

1 3 

Poo = -(ejj - -ej-ej) 

7 

t = e(4eu - -eid) 

For the operator acting on A we obtain: 
9ij = e<%, goo = 1 
r^- = ^{eiS jk + ej5 ik - e k 5ij) 

uf b = l;{-e a 5 b i + e b 5 a i - ^ei5 ab ) 

1 5 

E ab = M ab + — (e kk 5 ab + e a e b + -e p e p 5 ab ) (24) 

R l jki = -{ejAk — ej k 5u — eu5 k j + ei k 5ij) 

+^(e p e p (5ji5 ki - 5j k 5u) + e k ejdu - e k ei5ji - e\ejb ki + eierfjk) 

Pjk = ~nifijk Cppfijk) ~^(.^k&j CpCp^fcj) 

1 1 

r — — {— 2e pp — —e p e p ) 

Here for convenience we prefere to keep the distinction between coordinate 
indices k, 1} and bundle indices {a, b}, though they all run from 1 to 
3. In the equations (22) - ( |2~4| ) repeated indices are contracted with the 



flat space metric Sij. 

It is instructive to express the heat kernel coefficients in terms of e and 
its derivatives: 

*W/.t) = J 

t 2 , „„ „ 237 



_ t _ 2gg( 336jjj'j 186j6}j'j 336jj'6jj + ^ ^H^jj 

531 33 837 3 

H q Cij&idj H ^^ii&j&j ^ g^~^*^*^i^i) ~l~ C/(i J) 

13 

t 2 e 2 119 
~i~ ( 276jjjj 60ej6jjj 41ejjCjj H ~^~&a^jj 



91 , 415 4141 3 

^(/,t) = (^/^/{i+Ki^-^) (25) 

t 2 , _ 711 

360e 2 llleieijj + lo2eijeij e a e jj 

1029 793 | 4263 3 

^ dijGiCj H -^-CiiCjCj H — 6i€i€jCj) -\- CJ\t )} 

Here ej...j = di . . . djlne. This completes the calculation of the UV divergent 
terms. 

We can define a "total" heat kernel as i^[A ] + ^pji ~~ 2K g h- We see, 
that the contribution of ghosts is not cancelled by that of Aq and of the 
"non-physical" components of A. 

As a check, in the Appendix we derive (^) by an alternative method. 

The asymptotic expansion constructed above gives 2n spatial deriva- 
tives of e in any a n . Hence it is clear that certain smoothness of e(x) is 
needed. Our expansion is not valid if e changes abruptly, as e.g. for a 
bubble in water. For the configurations of latter type boundary terms in 
the heat kernel expansion should be taken into account. 



4 Conclusions and discussion 

In the present paper we performed the path integral quantization of elec- 
tromagnetic fields in a dielectric medium. As a first step, we considered 
the first order action and derived the canonical Poisson brackets. Next, 
we constructed the canonical (simplectic) measure in the phase space. We 
built up a measure in the configuration space by means of an integration 
over the canonical momenta. This measure appeared to be different from 
the naive one. By choosing a suitable gauge fixing condition (0) we re- 
duced the path integral to a product of three determinants of operators 
of Laplace type. For the evaluation of the ultra violet divergent parts of 
this determinants the standard heat kernel technique |l3| is available. Our 
results are re-checked by another technique (see Appendix). We observed 
no cancellation of ultra violet divergencies between ghosts and any "non- 
physical" components of the vector potential. Thus it is highly unlikely 
that the full quantized electrodynamics in dielectric media is equivalent to 
a theory where only two polarizations of photons are quantized. 

The next step to do is to work out a suitable cut-off procedure for the 
path integral. This problem is very non-trivial in the present case. Since 
e — * 1 at high frequencies, the cut-off is physical, it will not be removed after 
a renormalization. Therefore, we must be sure that the basic properties of 
the quantum field theory, as unitarity and absence of gauge anomaly, are 
valid at finite cut-off. After having solved this problem, it will be possible 
to consider the vacuuum energy densities and other physical quantities of 
interest. 
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Trexp(-tL) = J < x\ exp(-£L)| k >< k\x > (26) 



Appendix 

In this Appendix we describe briefly an alternative method for the evalu- 
ation of the heat kernel expansion which we used to control our results. 

We can represent the functional trace in the r.h.s. of ( |l~5|) as an integral 
over x of diagonal matrix elements between < x\ and \x > and insert 
"unity" expressed via an integral of momentum eigenstates: 

d 4 xd 4 k 

w? 

The generic form of the matrix element in ( p6| ) is < x\F\ (e i de)F2(d)\k >, 
where F\ and F2 are some polynomials of e and its derivatives and of 9$ 
respectively. Acting on the left F\ is replaced by its value in the point x. 
Acting on the right, F2 is replaced by F2(ik). It is easy to see that the 
result is 44 

/ exp(-iL(e(x), d M -> <9 M + (27) 

where we should take all external fields in the point x, shift all derivatives 
by ik, and drive derivatives to the right. It is understood, that d standing 
at the very right position vanishes. 

Consider the heat kernel for the ghost operator: 

K gh (t) = J ^0^- exp(t($ + 2ik j d j + 

+(dj log e)dj + ikj {dj log e) - k 2 - ecu 2 )) (28) 

where {k^} = {uj,kj}. Time derivatives are dropped out because e(x) is 
static. 



To obtain a small t asymptotic expansion of (28), one should isolate 
the factor exp(— t(k 2 + euJ 2 )) and expand the rest of the expression in a 
power series of operators and functions involved. Next one should integrate 
over momenta and collect all terms with the same powers of proper time t. 
Denote the exponential in ( [28] ) as exp(^l + B) , where A = —t(k 2 + ecu 2 ). 
Note, that A does not commute with B. However, the repeated commu- 
tator [[[B, A], A], A] vanishes. This allows us to present the exponential as 
follows (see e.g. Q) 

exp(A + B) = expA(l + B + l[B,A] + ^[[B,A],A]+ (29) 

+ \B 2 + \[B, A]B + l[B, [B, A]] + ^[B,A] 2 + ... 
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We retained all the terms which contribute to the two leading terms of the 
asymptotic expansion proportional to t~ 2 and f 1 . 

Acting as explained above we obtain the asymptotic expansions for the 
heat kernels K g h , ^[a ] an d K[Ay ^he nrs ^ ^ wo terms are m complete 



agreement with (25). Calculations of the third terms are too complicated 



to be done just for a control. 
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